Shortest paths on systems with power-law distributed long-range connections 
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We discuss shortest-path lengths i(r) on periodic rings of size L supplemented with an average 
of pL randomly located long-range links whose lengths are distributed according to Pi ~ l"'^. 
Using rescaling arguments and numerical simulation on systems of up to 10^ sites, we show that 
a characteristic length ^ exists such that i{r) ^ r for r < ^ but l{r) ^ for r » ^. For 

small p we find that the shortest-path length satisfies the scaling relation l{r, fj,,p)/£, = /(MjJ'/C)- 
Three regions with different asymptotic behaviors are found, respectively: a.) fj, > 2 where 9a ~ 1, 
b) 1 < /i < 2 where < Os{fi) < 1/2 and, c) ^ < 1 where i{r) behaves logarithmically, i.e. 9s = 0. 
The characteristic length ^ is of the form ^ ~ p"" with i/ = l/{2 — fi) in region b), but depends on 
L as well in region c). A directed model of shortest-paths is solved and compared with numerical 
results. 

PACS numbers: 05.10.-a, 05.40.-a, 05.50.-fq, 87.18. Sn 



INTRODUCTION 



It has been known for long that slowly decaying long- 
ranged (LR) interactions can drastically change the crit- 
ical behavior of a system. A well studied example 
is the onc-dimcnsional Ising model with J(r) ^ r^^ 



il. 

Kondo problem llll, 



^ which is relevant for the 
121 among others. If /i > 2 there is 
no ordered phase at any finite temperature, the same as if 
only short-ranged interactions were present. When /i — 2 
the magnetization undergoes a finite jump at > 0, 
while all derivatives of the free energy remain finite (es- 
sential singularity). When /i < 2 the model displays a 
second-order phase transition with /i-dependent critical 
indices, which take their classical, or Mean-Field (MF) 
values for ^ < 1.5. On approach to = 2 from below, 
the correlation- length exponent diverges, signaling the 
appearance of an essential singularity. This divergence is 
of the form ~ (2 — /i,)~^/^ Q for Ising and (2 — /i)~^ 
for n-component models with n> 1 (but see p3| , where 
V ~ [2 — ix)~^\/n is suggested). A comprehensive account 
of what is known for Ising systems with LR interactions 
has been given by Luijten and Blote ||]. 
For d-dimensional n-component systems with ferromag- 
netic interactions decaying as l/r''+°'. Fisher, Ma and 
Nickel in propose that the lower critical decay rate is 
given by (T = d/2, or equivalently that the upper crit- 
ical dimension is du = 2cr. For a < d/2 the critical 
indices take their MF values, for (i/2 < cr < 2 they are 
cr-dependent, and for a > 2 they take their short-range 
(SR) values. Similar investigations have been conducted 



for Potts II 0, H 111, HeisenberglH |T|, |l[ H, 
p^ , and other p4 I^Sj models. 

The following picture is often found: for small enough 
decay rate /x, MF indices are obtained. Upon increasing 
/i a regime follows where critical indices change contin- 
uously with /X until finally SR indices are recovered. In 
a loose sense one can say that the addition of LR in- 
teractions changes the "effective dimension" of the sys- 
tem, although in a way that may depend on the specific 
model considered. This idea has been exploited to study 
the scaling behavior of critical systems above their ef- 
fective upper critical dimension du, while still working 
on lattices of low Euclidean dimension . The connec- 
tion between LR interactions and dimensionality was also 
briefly touched upon by Scalettar . A possible way to 
define an effective dimension, which is in general model- 
dependent, is to do so through the hyperscaling relation 
(2 - a) = dv, as pi [l7| 4// = v-^{2 



10], [17] deff = V ^{2- a] 
An alternative paradigm for the problem of LR interac- 
tions considers systems on a d-dimensional lattice sup- 
plemented with randomly distributed LR bonds of unit 
strength, which are present with probability pij ~ '"'^j^ ■ 
Notice that in this case the system has disorder: it is the 
probability for a given bond to be present, and not its 
strength, what decays with distance. These two ways to 
introduce LR interactions; decaying strength (DS) and 
decaying probability (DP), are not in principle equiva- 
lent. It is well known that disorder may change the crit- 
ical behavior if the specific-heat exponent a is negative. 
The DP paradigm is on the other hand relevant for a 
number of problems in which connectivity, and not the 
strength of the interaction, is determinant of the physi- 
cal behavior. Examples of problems of this kind are the 
magnetic |2^ and conductive ^ properties of poly- 
meric chains, where the probability of crosslinks between 
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two monomers decays as a power-law of the chemical dis- 
tance between them, conduction in insulating matrices 
with one dimensional conducting inclusions [ pOl w hose 
length distribution is "broad", neural networks |3l|, 
geodesic propagation on spaces with topological singular- 
ities (wormholes), the spread of fire or diseases 
etc. 

Networks built according to the DP paradigm of LR in- 
teractions may be characterized entirely in geometrical 
(or topological) terms, because all bonds have the same 
strength. Thus it appears for example possible to define 
the relationship between effective dimension dg/ / and de- 
cay rate /i of interactions in purely geometric terms for 
these networks. 

A useful topological characterization of random networks 
is the Graph Dimension dg, defined as follows: if V{£) is 
the average number of sites that can be reached from a 
given one in £ steps between connected neighbors, then 
V{i) ~ f'^f asymptotically. We now let £{r) be the av- 
erage smallest number of links needed to join two points 
separated by an Euclidean distance r (the "shortest-path 
length"), which behaves asymptotically as where 6s 
is the shortest-path dimension Since V{r) r'', the 
above relations imply that dg = d/9s, and we see that the 
asymptotic behavior of i{r) defines the graph dimension 
dg. 

In this work we study shortest-paths on DP networks, 
i.e. d-dimensional lattices with the addition of an aver- 
age of p LR bonds (or shortcuts) per site, whose length 
is distributed according to Pi ^ l^'^. We shall con- 
centrate mostly on the case d = 1, where numerical 
simulations are easiest. DP networks with power-law 
distributed LR bonds have been recently considered in 
one dimension both from the point of view of Random 
Walk properties and Shortest-Path lengths [Q, but 
for small system sizes. We will later discuss some of the 
conclusions in [^ , which appear to need revision in the 
light of our results. 

In Section II several definitions which are relevant for our 
problem of shortest-paths on Id DP networks ar e giv en. 
Simple rescaling arguments are used in Section II A to 
show that /i = 2(i is a critical decay rate, such that for 
fi > 2d, LR bonds are unimportant on large scales. For 
fi < 2d on the other hand, when p is small these ar- 
guments predict the existence of a characteristic length 
£, ~ ^ beyond which LR bonds are important. In 

Section |IIB| a directed model is introduced for shortest- 
paths in Id, which turns out to be exact for /i > 2 and 
still provides an useful upper bound when /i < 2. In Sec- 



tion III our extensive numerical results for shortest-path 
lengths £{r) in one dimension are described and compared 
to theoretical predictions. Finally, Section IV contains a 
discussion of our results. 



II. DP NETWORKS AND RESCALING 

We start with an arbitrary c?-dimensional lattice made 
up of = L'^ sites, and its corresponding SR bonds. In 
addition to these, DP networks are defined to have an av- 
erage of pL'^ LR bonds, or shortcuts, whose lengths and 
locations are random. This is done in practice by let- 
ting one LR bond stem from each site i with probability 
p. The neighbor j at the other end of each LR bond is 
randomly chosen with a probability P{j\i), that is a de- 
caying function of the Euclidean distance = \xi — Xj\ 
between sites i and j. 

For a given realization of shortcuts, the shortest- 
path length £ij is defined as the minimum number of 
connected- neighbor steps needed to join sites i and j. 
This quantity is measured as a function of Euclidean 
distance r^ , and averaged over disorder (realizations of 
shortcuts). After disorder average, £{r) is the average 
"cost" of joining two points separated by an Euclidean 
distance r, and is defined as 

£{r) = ^ < % > 5{r,, ~ r)/ ^ S{r,, - r), (1) 



where <> means disorder average. 



A. Rescaling 

Consider now dividing the d-dimensional lattice into 
"blocks" of linear dimension b, such that 1 << 6 << 
L, and regard two sites I and J of this new lattice to 
be connected by a shortcut if any pair {i G I,j G J} 
is connected by a shortcut. We allow for at most one 
shortcut between rescaled sites since, for the purpose of 
shortest-paths, the only fact that matters is whether two 
sites are connected or not. If the original pairs ij are 
connected by a shortcut with probability pij , the rescaled 
probability qij = 1 — pu for blocks I and J not to be 
connected is given by 



QiJ= n ('^-Pv)= n 



(2) 



ieijeJ 



ieijeJ 



which for large distances » b can be approximated 

as qij = qlj ■ This can be written as q{r/b) — {q{r)) 
and therefore X{r) = log<7(r) transforms in a simple way 
under rescaling. 



X{r/b) = b'^X{r) 



Thus 



p{r) 



-p/r'- 



(3) 



(4) 



retains its functional form under rescaling, i.e. 



p(r) = (l - e-'^/'-^) , 



(5) 
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with p = h'^'^'^ ^^p. The condition that the system con- 
tains a total of pU^ LR bonds is ensured by imposing 



P - 



(6) 



where Sd is the surface of a c?-dimensional hypersphere 
of radius one. This relationship fixes p as a function of 
p and L. In the limit of small p, p turns out to be pro- 
portional to p. Notice that, because of the multiplicative 
rescaling Eq. a pure power law is not strictly invari- 
ant under rescaling. But the true invariant distribution 
Eq. can be very well approximated by a power law 
for large distances r such that p/r'^ << 1. Restricting 
ourselves to the limit of small p (or p) we can thus work 
with a power-law distribution of shortcut lengths. In the 
following we consider 



(7) 



where the normalization constant C is chosen so as to 
satisfy Eq. (||) . In Appendix A2 we show that p rescales 
as 



with 



p — b^pp, 



d for p < d 

2d — p for p > d 



(8) 



(9) 



Notice that expressions similar to Eq. and Eq. (||) 
give the renormalized coupling constant of the one- 
dimensional LR Ising model at low temperatures 
It follows that p = is a line of fixed points in the p,p 
space of parameters. For p < 2d this fixed line is repul- 
sive, and becomes attractive for p > 2d. Thus for p > 2d 
the density of LR bonds is renormalized to zero under 
rescaling, and Pc = 2d is the upper critical decay rate 
above which LR bonds are irrelevant, and SR behavior 
is recovered. 



B. Naive Paths: An approximate model in one 
dimension 

Consider a directed path which starts at t — from 
xq — 0, proceeds always to the right, and is built by 
using at each site any LR bond available, provided this 
bond does not take the path further to the right than 
r. We call the path so defined the "naive path" between 
and r. As compared with the actual shortest path, 
this construction neglects the possibility of turnbacks, or 
that certain LR bonds may not be used (See Fig. |^). 
We will later see that under certain circumstances, the 
naive-path approximation gives a reasonable estimate for 
shortest-path lengths. But even if this is not the case, 
the former constitutes an upper bound for the shortest- 
path length, and thus still provides useful information. 
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FIG. 1: Full lines are SR bonds, dashed LR bonds. One 
possible shortest path between and 8 is {0-1,1-5,5-4,4-8} 
and has length four. The naive path uses all rightwards LR 
bonds available at each site, i.e. {0-2, 2-3, 3-5, 5-6, 6-8} and 
has length five in this example 



The naive-path length is the number of timesteps 

it takes to reach r, and can be estimated in the following 
way. At time t the walker sits at site Xt- From this 
site, with probability p a LR bond (of random length l^) 
stems rightwards. The walker now proceeds along this 
LR bond, provided it does not go further to the right 
than r. The joint probability pt that a bond is present 
at Xt , and its length is not larger than r — xt is 



i — Xt 

Pt=P^ Pi 



(10) 



Thus at time t the walker goes one unit to the right with 
probability qt = 1 — pt, and It units with probability pt- 
Therefore in average 



Xt = Xt-l + 1 +pt{k - 1), 



(11) 



where It is the average length of a LR bond which is not 
larger than r — xt, i.e. 



p 



1=1 



Thus Eq. (y) reads 

Xt = xt-i + 1 + pG{r - Xt), 

where 



(12) 



(13) 



(14) 



Within a continuous-time, continuous-space approxima- 
tion we put 



x{t) = l+pG{r-x{t)) 



(15) 



which shall be solved with boundary conditions x{t — 
0) ~ and a;(t = £„(r) — 1) = r — 1 (notice that Eqns. 
( |l0| ) and ( p^ are only defined for Xt < r — 1). This can 
be formally integrated to give 



Lir) = 1 + 



dx 



1 l+pG{x)- 



(16) 



We will analyze this result and compare it with our nu- 
merical results in the following sections. 
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III. NUMERICAL RESULTS IN ONE 
DIMENSION 

In this section, numerical results are presented for pe- 
riodic rings of up to 10^ sites. One LR bond stems from 
each site with probability p < I. Its random length I 
is obtained by first generating a real random variable 
z such that 1 < z < {L/2 + 1) with P{z) - z~^, 
and then taking its integer part: I — Int{z). Lat- 
tice sizes are Lfe = lO^+'^Z^ foj. ^ 0,l,---,8. The 
density of LR bonds is p = 0.001,0.003,0.01,0.033,0.1 
and 1.0. Shortest paths are identified by Breadth-First- 
Search (BFS) ||3^, |39|], and averages are taken over 10^ 
samples. Altogether the results presented in this work 
involve an amount of computational work equivalent to 
approximately 10^^ sites. Figures I, I and I show av- 
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FIG. 3: Same as Fig. 0, for 1 < /i< 2. 



separated by an Euclidean distance r to be connected by 
a LR bond does not depend on L when ^ > d. (See 
Eq. (|A6D). In comparison, when fi < d one has that P{r) 
decays to zero with system size as L"'''^''^. This scale- 
dependence in the connectivity properties is evidenced 
by the size-dependence of i{r) when /i < 1 in Fig. ||. 
A second noticeable feature is that for all /i < 2 a char- 
acteristic size ^ exists with the following property: For 
r << ^, i{r) w r, while for r > ^, i(r) grows asymptoti- 
cally slower than r; in general as r^" with 9s < 1. This 
characteristic size ^ is a function of p and for 1 < /i < 2, 
but also depends on L for < 1. 



FIG. 2: Average shortest-path length £{r) vs r. Nu- 
merical averages (full lines) over 10* samples are shown for 
systems of size Lfc = iQ3+k/2 ^^^^i k = 0,1, •••,8. The 
dashed line is £{r) = r. The local density p of LR bonds 
is p = 10""^, 10~^, 10~^ and 1. The different cases can be told 
apart by noticing that larger values of p result in lower values 
oie. 



erage shortest-path lengths £{r), respectively for the re- 
gions: 0</^<l, l</i<2 and n>2. 
It is apparent in these plots that £{r) does not depend on 
system size L (only on p and fi) for /i > 1 . This is consis- 
tent with the fact that the probability P{r) for two sites 



A. The p. > 2 regime 

As seen in Section II, for /i > 2 the density of LR 
bonds rescales to zero, i.e. p = is an attractive fixed 
line. Thus one does not expect LR bonds to modify the 
effective geometry of the lattice in this regime. In fact 
it is found (Fig. ^) that £{r) oc r at large distances, and 
thus deff = d in this regime, although the coefficient 
of proportionality depends on /i and p in general. Our 



directed model (naive paths) described in Section [IB 
gives exact results in this regime as we now show. 
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FIG. 4: Same as Fig. |, for ^ > 2. 



1. Naive paths when /i > 2 



When fi > 2, G{x) in Eq. grows monotonically 
from G(l) = to G{oo) = 1 — 1. Thus asymptotically 
(-nir) = + — 1)). In order to obtain the short- 
distance behavior we may approximate, to first order in 



[l+pG{x)]-' ^l~pG{x) 
Eq. ( |l6| ) now reads 

^„(r) «r[l-p$(r)], 

where 

$(r) = - [ G{x)dx 



(17) 
(18) 

(19) 



is a p-independent function which converges to I — 1 for 
large r. Equation can be integrated (See Section |A^ ), 
and the comparison between analytical and numerical re- 
sults is done in Fig. ^ The coincidence betweeen the 
naive-path model and numerical results is very good even 
at /i = 2. Thus we conclude that in the fi > 2 regime 
and when p is small, shortest-paths are essentially naive 
paths. 



B. 



The 1 < fi < 2 regime 



In Section II A we saw that p 



is a repulsive 
fixed point for all /i < 2 in one dimension. Because 
of the rescaling law Eq. (||), one expects a lengthscale 
f ~ p~^lyp = ]-,g relevant for the behavior of 

£{r) as p — > 0. For r << ^, the p = fixed point is dom- 
inant (for which £{r) = r) while for r >> ^ the effects of 
LR bonds may become visible {i{r) shorter than r). 




FIG. 5: For fi > 2 shortest path lengths ^(r) are well approx- 
imated by Eq. ( |l8| ) with $ given by Eq. (^. Shown in this 
plot are our numerical results (solid lines) for $(r) = p~^(l — 
e{r)/r) for small densities of LR bonds: p = 10"^, 3 x 10"^ 
and 10~^ . The dashed line indicates our analytic result. 



Eq. (AlC) 



1. Naive paths when 1 < /i < 2 

For fi < 2, I is not well defined. However the average 
length /( of a LR bond not larger than r — xt is well 
defined and given by Eq. (p^). Notice that G{r) now 
grows as r^~^. Eq. ([l6| ) is still valid for naive paths, 
and one gets in the limit of large r that £{r) ~ r^~^, i.e. 
gnaive _ ^ _ turns out that actual shortest-paths 

are shorter than naive paths for /i < 2, i.e 6*"°™'^ = /i — 1 
is only an upper bound for dg (see Section III B 3| ) . 

Although the naive-path model fails to predict the 
asymptotic behavior of £{r), it can nevertheless still help 
us determine the characteristic length ^ beyond which 
£{r)/r 0. Keeping just the fastest-growing term 
in G{y) (Eq. (Ill)) and equating pG{^) ~ 1, one gets 



-1/(2-m) 



in full accordance with rescaling argu- 



ments in Section II A and at the beginning of this section. 



We show next that this is verified numerically. 



2. A single characteristic length £, 

In this subsection we test the hypothesis that a single 
lengthscale £,{p) dictates the behavior of £(r) in the limit 
of small p, and show that for 1 < fi < 2 this lengthscale 
is £, in accordance with rescaling arguments 

(Eq. (||)) and naive-path predictions. We propose that, 
for p ^ 0, 



^(^,M,P)/e = /(M,VO, 



(20) 
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where ^ p , and 

X for a: << 1 

/(m,x)oc<| (21) 
2,9. (m) for x»l 

This means that all p-dependence of l{r) is contained in 

By comparison with our numerical results we find 
that f{x) can be well approximated by f{x) = x/[l + 
Ca;(i-'^=)1. Therefore 



£(r)/e 



l + C (r/Oi- 



or, equivalently 



£{r) 



(22) 



(23) 



provide a good approximation to our numerical results. 
We fit Eq. ( p3[ ) to our numerical data for L = 10'' and 
p = 0.001,0.003,0.010 simultaneously (using v{fi),es{fi) 
and C(/i) as fitting parameters), and find and 6s as 
shown in Fig. |[ These results are entirely consistent 
with l/i/ = (2 — /i) for 1 <u < 2. Larger values of p are 
found not to follow Eq. ( po| ) satisfactorily, therefore we 
must regard this scaling expression as only valid in the 
p limit. 

A plot of ({r)/£,{p) vs. r/^ is shown in Fig. ^ for p = 
0.001,0.003 and 0.010. The fact that all three values of 
p collapse neatly onto one single curve suffices to verify 
the correctness of our scaling ansatz Eq. (pcl) f or small 
p. The specific form of f{x) chosen in Eq. (|22D should 
however only be regarded as empiric. 

Although for /i < 1 we do not expect Eq. (|2^) to 
hold (since then ^ has an additional L-d epend ence not 
included in these expressions, see Section IIIC ), a fit of 
the data gives ^ « 1, indicating that the p-dependence of 
the characteristic size ^ is of the form p~^ in this region. 
This is again consistent with Eq. (||) . We wil l discuss the 
regime n < 1 in detail later in Section III C . When p is 
small and fj, is close to 2, ^ grows too large. Consequently 
neither ^ nor can be correctly estimated for fi > 1.6. 
Consider for example p — 10~^. One then has ^ ~ 10^ 
for ^ — 1.6, but ^ ^ 10^*^, well beyond our present reach, 
for ^ = 1.8. Thus the estimates for 6 and 1/v in Fig. ^ 
are to be disregarded for /x > 1.6. 



3. Asymptotic exponent 9s 

When r >> we find that t{r) grows asymptotically 
as r®= . The shortest-path dimension 9s depends on /i 
only, goes to zero as /i ^ 1+ and jumps discontinuously 
to 9s = 1 sX ^ — 2" . We estimate 9s by two different 
methods. A simple power-law fit of the large-r behavior 
of £{r,fj,,p) gives the estimates shown in Fig. || for L 
ranging from 10^ to 10^ and several values of p. Strong 




0.4 0.6 0.8 1 1.2 
Decay Exponent 

FIG. 6; Numerical estimates for l/v (asterisks) and 9 
(plusses), obtained by fitting Eq. ( ]2^ ) to our data for L = 
lO'^ and p = 0.001, 0.003 and 0.010. The dotted lines are 
l/v — 2 — fi and 1/^ = 1. 




FIG. 7: Data collapse of l{r), showing plots of l/^ vs r/^ 
with ^(p) = p-^/f^-^'', for p = 0.001,0.003 and 0.010. The 
dashed line is our approximate expression Eq. (p^. For larger 
values of /i < 2 and the same values of p, the characteristic 
size is much larger than 10^. 



finite-size corrections affect the smaller values of p, for 
which >> L when fi ^ 2. However for large L all these 
estimates are seen to converge to similar values within 
numerical accuracy. 

The second method chosen to estimate 9s consists 
in fitting our numerical data using Eq. ( p3| ) but with 
V = 1/(2 — /i) instead of taking v as a fitting param- 
eter as in Fig. ^. Fits of our data for L ~ 10^ and 
p — 0.001,0.003,0.010 produce the values of 9s shown 
in Fig. ^. Again the results obtained for fi > 1.6 are to 
be disregarded since ^ is much larger than L for these 
values of p. A naive interpretation of the results in 
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FIG. 8: Asymptotic exponent 6s obtained from power-law 
fit of tiie large-r behavior of (.{r), for L of the form Lk ~ 
■|^Q3+fc/2^ fc = 0, 1 . . . , 8 , and for the values of densities of LR 
bonds p indicated in the respective plots. Lines are guides to 
the eye. 



Fig. g, for any fixed value of L, could lead one to believe 
that the transition between linear behavior {£{r) oc r) 
and sublinear behavior (£{r)/r — > for r oo) hap- 
pens at a p-dependent boundary ndp) |^7| . However, a 
more careful numerical analysis shows that this transi- 
tion happens at /ic = 2 for all p in the thermodyn amic 
limit, as predicted by rescaling argu men ts (Section p A| ) 
and the naive-path model (Section |ll B|) . This appears 



0.8 
0.6 
0.4 
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1 1.2 1.4 1.6 1.8 2 

Decay Exponent 

FIG. 9: Numerical estimates for the asymptotic exponent 
6s (plusses) in Eq. (p2[), resulting from fits of our data for 
L = 10'^ and p = 0.001, 0.003 and 0.010 with ^ = p"^/'^""'. 
The rightmost two points, for ^ larger than 1.6, suffer from 
strong finite-size effects and should be disregarded. The dot- 
ted line sketches what we believe is the true value of 6s{i-l). 
The discontinuity at /i = 2 is suggested by the behavior of 
the p = 1.0 results in Fig. ^. 

to be in partial disagreement with recent work of Sen 
and Chakrabarti (SC) , where the "regular lattice be- 
havior" {£{L) ~ L) is claimed to extend below ^ — 2 



for small values of p. SC explain what they call the 
lack of small-world behavior in lattice polimers as being 
a consequence of the small number of LR connections 
(small p). Based on the analysis of £{L) on relatively 
small (L — 10^) systems, SC conclude that there is a 
p-dependent phase boundary fidp) < 2, and show that 
several lattice polymer models lay marginally on the reg- 
ular lattice {£{L) oc L) side of this boundary. Our exten- 
sive numerical results and analytic considerations how- 
ever show that ^ — 2 is the critical decay rate below 
which £(r) << r, for any density p of LR bonds. The 
p-dependent boundary that SC observe is just a logarith- 
mically slow finite-size effect. At sufficiently low values 
of p, and for /i close to but lower than two, the charac- 
teristic length ^{HjP) is larger than L and thus £{L) cx L. 
Equating ^ = = one obtains an apparent 

boundary p*{p) = 2 — log(l/p)/ log(i), which converges 
logarithmically slow to pc — 2. Replacing L = W^, this 
last expression follows closely the boundary reported by 
SC in Figure 3 of [0. 

There is a second aspect of with which our find- 
ings seem to be in disagreement. According to SC, there 
are only two phases regarding the asymptotic behav- 
ior of £{r). A logarithmic phase, £{L) cx log(L), for 
/i < /i*(p) ~ 2, and a linear phase for /i > fJ.*{p). Our nu- 
merical evidence however suggests a more complex sce- 
nario. For 1 < /i < 2 we find that £{r) oc with dg 
small but nonzero (Fig. |£ ), an d only for p < 1 £ becomes 
logarithmic (See Section IIIC). 



C. The < fi <1 regime 

The data in Fig. |^ clearly show that £{r) depends 
on system size L if < /i < 1. In the specific case 
/I = 0, each of the L^'^/2 possible LR bonds is present 
with the same probability pL^"^. This corresponds to a 
d-dimensional lattice supplemented with pL''' LR bonds 
whose ends are randomly chosen, and goes under the 
name of Small- World (SW) network ||, ||, 
H, ||,J4^, E6[ 1^, Hi- In particular it was recently 
found |33| , that on SW networks (/i — 0) there is 
still a single characteristic length Vc dictating the behav- 
ior of £{r), but it depends both on L and p, and diverges 
as L — > cxD, in any dimension d. Analytic calculations 
confirmed by numerical measurement |^ , ^ show that, 
in d dimensions, 



£ir) 



r for r < Tr 



for r > 



(24) 



where Tc ~ p^^^'^log{KpL'^) with K a constant. 
In the particular case d — 1 one has rc{fi — 0,L,p) ^ 
p~^ log(4pL). So the fi — case is relatively simple, with 
£(r) = r for r < log(4pL)/p and £{r) = log(4pL)/p for 
large r. 

By inspection of Fig. one concludes that Tc depends on 
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10" 



10" 



|i= 00.000 



|i= 00.500 





FIG. 10: Shown is i'/rc vs r/rc, with 

= p"Mog(4(pL)(i-''))/log(4). The density of LR bonds 
is p = 0.001, 0.003 and 0.01. System sizes L are of the form 
Lfc = 103+*=/2 ^ith fc = 0, 1 . . . , 8. The proposed expression 
for Tc has only been justified theoretically for /i = [ ps] , ^ , 
and is purely empirical for < /j. < 1. 



L as well as on p, n in the whole < /i < 1 range. When 
/i = 1 however, the characteristic length dictating the 
behavior of i{r) is ^ = p~^, and no longer i-dependent, 



as shown in Section [II B. Guided by these observations. 



we now propose an empirical expression for rc in terms of 
L,p, p in the whole < /i < 1 range. This expression has 
to result in ^ = p~^ when — > 1, and Vc ~ p~^ \og{pL) 
when /z — > 0. It is easy to verify that 



r,{ii,p,L)=p-Hog f4(pi)(i-'')l /log(4) 



(25) 



satisfies both requirements. We find that this empiri- 
cal expression gives acceptable results for small p. In 
Fig. nOwe show £{r)/rc versus r/vc with rc given by 
(Eq. (p5|)), for all values of L ranging from 10'^ to 10^ 
and p = 0.001,0.003 and 0.010. The acceptable collapse 
of all data supports the validity of (Eq. (^)) reasonably 
well. 

We find that £{L) grows asymptotically as log(L) for 
< /i < 1. The naive path model already predicts a 
logarithmic behavior at /i = 1 as the following shows. 
For fi = 1 one has (see the beginning of Section III) 



=log((Z + l)/0/log(L), from which G(a-) « x/ log(i:). 
Thus Eq. (0) can be written approximately as 



dx 



1 + x/rc 



(26) 



where 



p ^log(L). Thus naive-paths are deter- 



mined, in the /i — > 1 limit, by a logarithmically L- 
dependent characteristic size and a logarithmic behav- 
ior l{r) ^ log(r) above Tc- Given that actual shortest- 
paths must be shorter than naive-paths, we conclude that 
l{r) is logarithmic for all /i < 1. 



IV. CONCLUSIONS 

We considered shortest paths on c?-dimensional lattices 
of L'^ sites supplemented with pL'^ long-range connections 
whose lengths I are random variables with power-law dis- 
tribution P{1) ~ We call these decaying probability 
(DP) networks, since it is the probability to have a LR 
bond of length I, and not its strength, what decays with 
distance. The limit /i — + is the "small-world" network 
of Watts and Strogatz |Q . Under a rescaling transfor- 
mation with scale parameter 6 in c? dimensions, a small 
local density p of LR bonds transforms as p = b^'^~^p. In 
the (/i,p) plane, p = is a repulsive fixed line for pL <2d 
and an attractive fixed line for /i > 2c?. Thus rescaling ar- 
guments predict /Zc = 2d to be a critical decay rate above 
which LR bonds are irrelevant. Particularizing to rf = 1, 
a directed model that gives an upper boun d for shortest- 
paths can be analytically solved (Section |ll E| ) and has 
three regions in the /i-axis: a) t{r) cx r for /i > 2, b) 
l{r) cx r'^~^ for 1 < < 2 and c) l{r) logarithmic for 
/i < 1. In accordance with rescaling arguments, we find 
numerically that in one dimension /i = 2 is a critical point 
separating a "short-range phase" (/i > 2) where shortest- 
path lengths are linear, i{r) oc r, from a "long-range 
phase" (/i < 2) where shortest-path lengths are sublin- 
ear, l{r) oc r^' with Og < 1. Our finding that fic — 2 for 
all p is consistent with previous work of Jespersen and 
Blumen but is in disagreement with recent claims 
of Sen and Chakrabarty who suggest the existence 
of a p-dependent boundary At*(p). We showed that this 
apparent boundary is a finite-size effect, due to the fast 
growth of a correlation length ^ as /j, — * 2^. 

For small p and 1 < /J. < 2, a characteristic size £, = p~'^ 
with — 1/(1 — ^) dictates the shortest path properties. 
For r < ^ one has £{r) « r while for r >> ^, i{r) ^ 

(/i) is found. This divergence in the correlation length 
exponent as ^ — > 2~ is of the same kind as reported 
for spin models previously |0, ||, 13 1. 

For fj, < 1 the characteristic size behaves as p~^ but 
is also L-dependent and we find that Eq. (p5| ) provides a 



good empirical fit of both its p- and i-dependence. 

The asymptotic exponent 6s is found numerically to 
attain its short-range value 9s ~ 1 for /i > 2. It is dis- 
continuous at /i = 2, where it probably takes a value 
near 1/2, and then goes to zero smoothly as /j, ^ l"*". 
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For /z < 1 we find logarithmic (or Mean Field) behav- 
ior: 9s — and £{r) ~ log(r') asymptotically. For ^ — > 
£(r) saturates at large distances to a value which depends 
logarithmically on system size Es, H . 
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APPENDIX A: SHORTCUT DISTRIBUTION 

1. Normalization 

The scale-invariant shortcut distribution p{r) Eq. (^) 
can be approximated by p{r) « 1 for r < Tc = p^^^ and 
p{r) « p/r^ for r > re- Thus the normalization condition 
Eq. (pf) can be written as 



P = Sd 



J^" r'^-^dr + p r'^-^-^'dr for > 1 



pJ^r'^-^-t'dr 



for < 1, 



(Al) 

so that ifVd = Sd/d is the volume of a unit radius sphere, 

(A2) 



P = Vdi 



iii^ - 1 for p > 1 



for p < 1 



When p < d and if p remains finite in the L ^ oo limit 
one has that 



p = p- 



d- P r-{d- 



Sd 



j^-{d-n) 



(A3) 



This goes to zero for large L, which justifies the power- 
law approximation Eq. (^. For p > d on the other hand, 
and assuming p small. 



p=p- 



p — d 
Sd 



(A4) 



so that the power-law approximation holds for any finite 
p when p < d but only for p small when p > d. The 
power-law distribution is properly normalized when 



1 = Cy r'^-^-^'dr 



{p ~ d)Lf^-^ 
^ Sd{L.-^ - 1) 



therefore in the limit of large L 
pir) = 



^ ^ iovp>d 



S^Ld-f- rf lor p < a 



(A6) 



gives the probability for two sites separated by an Eu- 
clidean distance r to be connected by a LR bond. 
2. Rescaling 

From the rescaling of p, Eq. (||), and the relationships 
( |A3| ) and (A4) between p and p it is immediate to con- 
clude that, 



ith 



Vp 



p = b^pp, 

d for p < d 
2d — p for p > d 



(A7) 



(A8) 



3. Naive paths when /i > 2 

For p > 1 and L » 1 we have that P; = — 
(I + 1)^^^. Using this expression, Eq. ( p^ gives G{x) = 
ELi (l ~ 1) [^'"^ - + =H{x^\,p-\)-\- 
(a;-hl)"^ + (a; + l)"(^"^\ where E(x,a) = YTi^x lA" are 
called Harmonic Numbers. ll{x, a) can be approximated 
(within one percent error) for all a > 1 and a; > 2 by 



H{x,a) « 1 + 2" 



(A9) 



Within this approximation one obtains (Z — 1) = 
i/(oo,/i- 1) - 1 = 2^"^ + + which is found to 
be very accurate for all /i > 2. Using this approximate 
expression, Eq. (|9h can be integrated to give 



$(r) = 



(r-l)(r-l) Z(r,3-p) 

r r(p — 2) 

3Z(r, 2 - p) _ Z{r, 1 - p) 
2r r 



(AlO) 



where Z{x, a) = {{x + 1)" - 2")/c 
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